Opposites Attract - A Theorem About The Casimir Force 
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We consider the Casimir interaction between (non-magnetic) dielectric bodies or conductors. Our 
main result is a proof that the Casimir force between two bodies related by reflection is always 
attractive, independent of the exact form of the bodies or dielectric properties. Apart from being 
a fundamental property of fields, the theorem and its corollaries also rule out a class of suggestions 
to obtain repulsive forces, such as the two hemisphere repulsion suggestion and its relatives. 



The Casimir effect has been a fundamental issue in 
quantum physics since its prediction Q. The effect has 
become increasingly approachable in recent years with 
the achievement ofprecise experimental measurements of 
the effect |a, Li Lll- LJ , probing the detailed dependence of 
the force on the properties of the materials, and measur- 
ing new variants such as corrugation effects. The theory 
and experiment have good agreement for simple geome- 
tries. 

In spite of the vast body of work on the subject (For a 
review see 6] ) , some properties of the force are yet under 
controversy. Due to the computational complexity of the 
problem, the main body of work on the effect is a collec- 
tion of explicit calculations for simple geometries. In this 
Letter we resolve one of these controversies and supply 
general statements about Casimir forces, applicable to a 
broad class of geometries. 

The interest in repulsive Casimir and Van Der Waals 
forces has grown substantially recently due to possi- 
ble practical importance in nano science, where such 
forces may play a role as a solution to stiction problems. 
It is known that repulsive forces are possible between 
molecules immersed in a medium whose properties are 
intermediate between the properties of two polarizable 
molecules 0- Conditions for repulsion between para- 
magnetic materials and dielectrics without recourse for 
an intermediate medium were given in jjj. However, the 
prospect of realizing materials with nontrivial permeabil- 
ity on a large enough frequency range is unclear [j| . 

It is common knowledge, based on the Casimir-Polder 
interaction, that small dielectric bodies interacting at 
large distance attract . Based on summation of two 
body forces one may speculate that any two dielectrics 
would attract at all distances. In this Letter we show that 
at least for the case of a symmetric configuration of two 
dielectrics or conductors this prediction holds indepen- 
dently of their distance and shape, for models which can 
be described by a local dielectric function. Of course, in 
any real material as distances become small enough, i.e. 
compared with interatomic distances, Casimir treatment 
of the problem is not adequate anymore. 

We first emphasize that the two-body picture is not 
enough to prove this. Calculations of the interaction be- 
tween macroscopic bodies by summation of pair - interac- 



tions are only justified within second order perturbation 
theory. Indeed, in [j| it was demonstrated how summing 
two body forces may give wrong prediction for the sign 
of interaction between extended bodies. 

Another objection to the pair-wise intuition is based 
on the example of Casimir energy of a perfectly conduct- 
ing and perfectly thin sphere. This was worked out by 
Boyer and yields an outward pressure on the sphere. 
This result motivated a class of suggestions for repulsive 
forces, the most well known of which are two conducting 
hemispheres - considered as a sphere split into two and 
therefore expected to repel each other 0,0 (Fig. QJ. 

One may try to use perturbative series, such as the 
multiple scattering series in the conducting case and 
show the attraction term by term. However, checking 
such a claim at orders higher then second might prove 
a difficult task. Such an approach is justified for distant 
bodies, but doesn't seem to be particularly promising for 
the problem at hand. 

Our main result is that the electromagnetic field (EM) 
or a scalar field, interacting with (non-magnetic) bodies, 
which are mirror image of each other and separated by 
a finite distance, will cause the bodies to attract. In 
particular, this shows that two hemispheres attract each 
other. The result holds for a scalar field in any dimension 
and even when the bodies are inside an infinite cylinder 
of arbitrary cross-section (perpendicular to the reflection 
plane) with arbitrary boundary conditions (b.c.) on the 
cylinder, thereby verifying and generalizing recent results 
for a Casimir piston 

Expressing the Casimir interaction as a (reg- 
ular) determinant. Several expressions are available 
for Casimir forces between dielectrics. We find the path 
integral method 0, 0, a convenient starting point 
for the presentation (alternatively, the result may be ob- 
tained using other approaches such as the Green's func- 
tion method) . We start with the case of a scalar field for 
simplicity, and explain later how the result is extended 
to the EM field. The action of a real massless scalar field 
in the presence of dielectrics can be written as 

S = \J^ d vJ ^C(V 2 + cA(x,a,))^ (1) 
where <j)* u — (j)-^, and e(u>, x) = 1 + x( x j w ) is the dielec- 
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FIG. 1: What is the direction of the force between two con- 
ducting hemispheres? While the outward pressure on a con- 
ducting shell might suggest repulsion, it follows from the ar- 
guments below that the hemispheres in fact attract 



trie function (we use units K = c = 1). The change in 
energy due to introduction of x m the system is formally: 

Ec — E x - E X=Q = 

-i / — logdet A (l + w 2 x(x,w)(V 2 + uj 2 + iOy 1 ) 
Jo 27T 

A determinant is mathematically well defined if it has the 
form det(l + A), where A is a "trace class" operator, i.e. 

| Aj | < oo with Xi eigenvalues of A (For properties see 
|2fj). The expression above is not of this form, and only 
has meaning when specifying cutoffs. Removing physical 
cutoffs will leave us with an ill defined determinant and so 
we keep in mind cutoffs at high momenta in the notation 
detA (one may use instead lattice regularization). 

At high frequencies x(o;,x) — > 0, provides a physical 



frequency cutoff. x( w ) an d (V 



- iO) are analytic for 



Rew, Imw > 0, justifying Wick-rotation of the integration 
to the imaginary axis iui ending up with: 



E ( 
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logdet A (l + w 2 x(x, i w)Go(x,x')) (2) 



Where G (x,x') = (x| 



|x'). Restricting the op- 



H A ®H B 

It turns out 



-V 2 +LU 2 

erator (1 + lu 2 xGq) to the support of x (more precisely 
to L 2 (Supp(x))) clearly does not affect its determinant. 
We assume x 1S nonzero only inside the volumes of the 
two dielectrics A, B and we therefore consider in the 
following (1 + uj 2 xG a ) as an operator on H A © H B — > 
H A ®H B where H A = L 2 {A) and H B = L 2 (B). It is 

then convenient to write it as (1 + lu 2 xGq 

1 A + oj 2 xaG oaa uj 2 xaGo A b 
u 2 xbGq BA Ib+u 2 xbGq BB 
that the part of the energy that depends on mutual posi- 
tion of the bodies, and as such is responsible for the force, 
is a well defined quantity, independent of the cutoffs. To 
see this, we subtract contributions which do not depend 
on relative positions of the bodies A, B: 



E c = E c (A\jB) E c (A) - E C (B) (3) 

As in this amounts to subtracting the diagonal con- 
tributions to the determinant which are not sensitive to 



FIG. 2: Bodies A and B are related by the reflection J 



the distance between the bodies, (i.e. only contributes to 
their self energies). This yields: 



E c = 



do; 
o 2^ 



| log det A 



1 + u 2 xaG 0aa uj 2 xaG 0ab 
u 2 XbG 0ba 1 + u 2 xbG 0bb 



, 1 + lo 2 XaG aa \ 1 



duj 



1 + lu xbG obb 
1 TaGqab 



(log det A . T n , 

ZTT \ 1 B LrQ BA I 



(4) 



where T a = Go — Finally, using the relation 
def ( } r * I = det(l — YX), which holds for block ma- 



Y 1 
trices we have: 

E c(o) = J °° % logdet(l - T A G 0AB T B G 0BA ). (5) 

Note that the (hermitian) operators T a are exactly the 
T operators appearing in the (Wick rotated) Lippmann- 
Schwinger equation [18|. Indeed, one may alternatively 
derive Eq. © within Green's function approach and us- 
ing T operators. 

In JSJ we disposed of the cutoff A as the expression is 
well defined in the continuum limit. We recall that an op- 
erator M is called positive (denoted M > 0) if (ip\M\ip) > 
for any tp. Since x(} w i x ) > ( as follows from general 
properties of the dielectric function 19]), the T opera- 
tors are positive and bounded, and T a Gq AB T b Gq BA is 
a trace class operator without need for cutoffs for any fi- 
nite bodies A, B [2^]. In fact, this holds also for nonlocal 
X as long as f(x) — > J A x(* w i x > x')f(x')dx' is a bounded 
positive operator H A — > H A . At this point the determi- 
nant is regularized and rigorously well defined for every 
lj and the integration over u> is convergent due to the 
exponential decay of the kernels Gq AB as lu — > oo. 

It is worthy to note that (for x > 0) all eigenvalues 
A of the (compact) operator T a Gq AB T b Gq BA appearing 
in © satisfy 1 > A > 

The Theorem. Having established a mathematically 
well defined expression for the Casimir energy, we now 
come to the main result: Consider a configuration of two 
bodies A, B related by a reflection (Fig. EJl, with x(*M) 
a bounded positive operator and separated by a finite dis- 
tance a; then (for fixed spatial orientations of the bodies) 
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Ec given in (JSJ is a monotonically increasing function of 
a (i.e. the Casimir force is attractive). 

Proof: We assume that A is located entirely in the 
negative x n half space, and that B is its mirror image 
under reflection through the x n — a/2 plane. To exploit 
the reflection symmetry we introduce a mapping J : A — > 
B given by J(x±,x n ) = (xj_, a— x n ). Note that B = J A. 
J is volume preserving and induces a unitary operator 
J : H A -> H B defined by Jip(x) = ip(J(x)). In the 
case that ip is a vector field, as in the EM case below 
we take Jip(x) = (^±(J(x)), -^ n (J(x))) (see Fig. 0. 
Since the bodies A, £? are related by reflection we have 
T B = JT A J f and thus: 



En = 



duj 
2^ 



logdet(l- T a G 0AB JT a J^G 0B a) (6) 



Note that GqabJ — J GqbA is a hermitian operator 
(this can be verified). The energy can therefore also be 
expressed as 



(7) 



Ec = / £ logd6t i 1 ~ (VTa~G oab JVt^) 



We now show by a direct calculation that (as operators 
on H A ): 



GqabJ > 
d a G 0AB J < 

Let 1(a) = (fp\Go AB J'\'ip} for a function ^(xj 
1(a) is explicitly given by 



(8) 
(9) 

,x n ) e Ha- 



1(a) = / dxdx' 

J Ax A 



dk 

(2^) 



^*(x)V(x') 



■ (x_l — )-\-ik n (x n -\-x f —a) 

x ^ s (10) 



k 2 



■ ur 



Note that x n + x' n — a < 0, allowing integration over k n 
by closing a contour from below the real k n axis: 



/(a) = / dxdx' 

JAxA 



dk. 



(2tt 



icf— 1 



dxi/'*(x)e i 
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(11) 



showing that /(a) > 0, which proves JSJ, and that 
d a I(a) < which proves ©. 

From (|8I9|I it immediately follows that the opera- 
tor Y = \/TaGqabJVTa '■ Ha — * Ha also satisfy 
Y" > 0, d a Y < 0. Hence a Feynman-Hellman argu- 
ment implies that all its eigenvalues 1 > A„(a) > 
are monotonically decreasing as functions of a. Since 



logdet(l — F 2 ) = 5Z„log(l — A 2 ) is absolutely conver- 
gent it follows c> a logdet(l — Y 2 ) > 0, and hence by Q 
also d a E c > 0. 

This completes the proof for the scalar case. 

To treat the EM case we start with the well known 
expression Eq. (80.8) of Lifshitz and Pitaevskii 0] for 
the change in free energy due to variation of the dielectric 
function e at a temperature T: 



5F = 5Fq + ]-T 



£ J>Tr{V5e). (12) 



Here Fo is the free energy due to material properties 
not related to long wavelength photon field, and w„ = 
2imT are Matsubara frequencies. T> is the temperature 
Green's function of the long wave photon field given by 

T>(x,x',tw)ij =< x\ VxVx+^ e (r,iM)l^' >ij 

EqlJBJ may be written as SF = SF + SF C where 



F c = h T J2Z-J l °Sdet A (W x V x +u;le(x,tLJ n )§l3) 

-logdet A (V x V x +u 2 n )] 
= l T Er=-oo lo S det A(l + w 2 x(a;,^„)X> (iw„)). 



where T> (x, x' , iuj n )i 



< x\ 



Note that 



)ij - ^IvxVx+uJ l x 

Fq is exactly the same as (0 , with the scalar propagator 
Go replaced by the vector propagator T> Q . 

Thus, starting with this expression, one repeats J3J and 
@ to get JSJ) , replacing Go by Vq everywhere (including 
in the definition of the T operators). The analysis of the 
determinant now proceeds exactly as in the scalar case. 
The only place in the proof which needs to be modified 
is where the explicit form of Go was used i.e. Eci. (|10fl . 
where we now have to use I?0y (k, iuj) — grxj2 (<% + -^-) 
instead. 

The effect of using the vectorial propagator in 
Eq.JUJ is to replace V*(x)^(x') by V 4 *(x)^(x')(% + 
^jf). In the vectorial case J acts by = 
(i/j±(J(x)), —ip n ( J(x))) so we get a factor (— l) s i n . Sub- 
stituting and integrating over k n as before we find 



fvec (^0 



dk 



(-1)^*^. + ^) 



(14) 



where <£ 3 -(kj_) = dxV> 3 (x) e - iki - x - L e a! "\ /k i+" 2 . Now it 
is straightforward to check that the expression in square 
brackets is positive for any <j>i and the theorem follows. 

Extensions and remarks: 
1) Finite temperatures: As remarked above we have 
/ Ijr — * at finite T. Since the positivity argu- 

ments apply to the determinant at each fixed imaginary 
frequency co, they will also hold at finite T. 
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2) Confined geometry in transverse directions Our the- 
orem is easily extended to cover the case when plac- 
ing the system inside an infinite cylinder, perpendicular 
to the x n — plane, with arbitrary cross section. In 
this case, one has to replace our Go by the appropriate 
Helmholtz green's function in the cylinder: Go(x,x') — 

J dk n J2j V3 ^+ll i +E J ± ) e ' fc " (x " } > where <Pn(x) are the 
appropriate quantized eigenmodes in the transverse di- 
rection, and the integration over kj_ is replaced by dis- 
crete summation. Substituting this expression in the rel- 
evant integrals such as (|1 If) yields the attraction. Since 
the attraction is independent of the ipj , this result is in- 
dependent of the b.c. one sets on the containing cylinder. 

3) Dielectric in front of mirror Suggestions were raised 
for re puls ion between arrays of dielectrics and a mirror 
plane [2JJ , based on results for a rectangular cavity. Vari- 
ation of our theorem shows that one actually has attrac- 
tion. Consider the body A to the left of a Dirichlet mirror 
located at x n = a/2. By the image method the propaga- 
tor is replaced by G(x, x') = Go(x— x') — Go(x+x' — ah). 
This may also be written as G — Go = —GqJ. It is 
then straightforward to arrive at the expression for the 
energv|2Rj| analogous to © with det(l — GqJTa) replac- 
ing det(l — GqTaGqTb)- Using similar considerations as 
in the proof above the attraction follows. 

4) Dirichlet b.c: Our approach never uses directly b.c. 
on the dielectrics; instead, we consider interaction with 
an arbitrary permittivity e{x,u>). This is adequate for 
describing real conductors. Idealized Dirichlet b.c. for 
a scalar field and ideal conductor b.c. for EM field are 
obtained as the limit of large \(i(J); However, Neuman 
b.c. do not follow from the present treatment, since they 
do not correspond to a positive perturbation, or indeed 
to any regular perturbation. 

5) Nonpositive perturbations. Cases of effective x < 
typically occur when the medium between the bodies has 
higher permittivity then the bodies. These cases as well 
as cases with nontrivial magnetic permeability \x may be 
covered in a way similar to the above theorem. How- 
ever conditions on x, /x must be specified to ensure that 
the eigenvalues of 1 — TaGqTbGq remain positive. These 
conditions are related to the assumption that the pertur- 
bation may not be so negative as to introduce negative 
energy modes into the system. 

Summary: Our main result is that the Casimir force 
between two dielectric objects, related by reflection, is 
attractive. Our theorem serves as a no-go statement for 
a class of suggestions for repulsive Casimir forces. Of 
course, the treatment is only valid at distances where the 
system may be described reliably in terms of the field and 
local dielectric functions alone. Although the above proof 
applies only to symmetric configurations, the approach 
presented here may be used to analyze the more general 
cases. A natural question rises: How far can our result 
be generalized? Which classes of interacting fields obey 



it? 
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